Velocity selection problem for combined motion of melting and solidification fronts 
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We discuss a free boundary problem for two moving solid-liquid interfaces that strongly interact 
via the diffusion field in the liquid layer between them. This problem arises in the context of liquid 
film migration (LFM) during the partial melting of solid alloys. In the LFM mechanism the system 
chooses a more efficient kinetic path which is controlled by diffusion in the liquid film, whereas the 
process with only one melting front would be controlled by the very slow diffusion in the mother 
solid phase. The relatively weak coherency strain energy is the effective driving force for LFM. As 
in the classical dendritic growth problems, also in this case an exact family of steady-state solutions 
with two parabolic fronts and an arbitrary velocity exists if capillary effects are neglected . We 
develop a velocity selection theory for this problem, including anisotropic surface tension effects. 
The strong diffusion interaction and coherency strain effects in the solid near the melting front lead 
to substantial changes compared to classical dendritic growth. 



The early observation of liquid film migration (LFM) 
were made during sintering in the presence of liquid phase 
Q or during partial melting of alloys (see ^ for a re- 
view). Nowadays LFM is a well established phenomenon 
of great practical importance. In LFM one crystal is 
melted and another one is solidified. Both solid-liquid 
interfaces move together with the same velocity. In the 
investigated alloys systems the migration velocity is of 
the order of 10~^ — 10~^ cm/s and it is controlled by the 
solute diffusion through a thin liquid layer between the 
two interfaces Q . The migration velocity is much smaller 
than the characteristic velocity of atomic kinetics at the 
interfaces. Therefore, both solids should at the interfaces 
be locally in thermodynamic equilibrium with the liquid 
phase. On the other hand, these local equilibrium states 
should be different for the two interfaces to provide the 
driving force for the process. It is by now well accepted 
(see, for example, 0i0|) that the difference of the equi- 
librium states at the melting and solidification fronts is 
due to the coherency strain energy, important only at the 
melting front because of the sharp concentration profile 
ahead the moving melting front (diffusion in the solid 
phase is very slow and the corresponding diffusion length 
is very small). Thus, the liquid composition at the melt- 
ing front, which depends on the coherency strain energy 
and on the curvature of the front, differs from the liquid 
composition at the unstressed and curved solidification 
front. This leads to the necessary gradient of the concen- 
tration across the liquid film and the process is controlled 
by the diffusion in the film. 

If only the melting front existed, the melting process 
would be controlled by the very slow diffusion in the 
mother solid phase and elastic effects would be irrelevant. 
In the LFM mechanism the system chooses a more effi- 
cient kinetic path which is controlled by the much faster 
diffusion in the liquid film. However, in this case the 
relatively weak coherency strain energy is involved as ef- 
fective driving force for this process. In this respect the 
LFM mechanism is similar to other well-known phenom- 
ena such as diffusion induced grain boundary migration 



and cellular precipitation In these processes a rela- 
tively fast diffusion along the grain boundaries controls 
the kinetics and the coherency strain energy plays also a 
controlling role. 

Thus, a theoretical description of LFM requires the 
solution of a free boundary problem for two combined 
moving solid-liquid interfaces with a liquid film in be- 
tween. In Ref. lal this problem was considered for simpli- 
fied boundary conditions: the temperature and the chem- 
ical composition along each interface were kept constant. 
Their values arc different for the melting and solidifica- 
tion fronts and differ from those far from the migrating 
liquid film. This means that any capillary, kinetic and 
crystallographic effects at the interfaces were neglected. 
It was found that under these simplified boundary con- 
ditions two co-focal parabolic fronts can move together 
with the same velocity. The situation is rather similar to 
a steady-state motion of one parabolic solidification front 
into a supercooled melt 0, Ql or one parabolic melting 
front into a superheated solid. In this approximation the 
Peclet numbers were found, but the steady-state velocity 
remained undetermined at that stage. Thus, the problem 
of velocity selection arises. 

Solvability theory has been very successful in predict- 
ing certain properties of pattern selecting in dendritic 
growth and a number of related phenomena (see, for ex- 
ample, 0,0,01). In the two-dimensional dendritic case, 
the basic approach is as follows. One attempts to model 
the dendritic tip by a needle crystal, that is, a shape- 
preserving steady-state growth shape which is a solution 
of the equation of motion governing diffusion in the neigh- 
borhood of a solidification front. This needle crystal is 
assumed to be close in shape to the parabolic Ivantsov 
solution. If anisotropic capillary effects are included, a 
single dynamically stable solution is found for any exter- 
nal growth conditions. This theory has been extended to 
the three-dimensional case 0, ^3 ■ We note that cap- 
illarity is a singular perturbation and the anisotropy of 
the surface energy is a prerequisite for the existence of 
the solution. 
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FIG. 1: Schematic presentation of two moving nearly 
parabolic fronts; 5*1 and S2 are the melting and growing 
solids, and L is the liquid film. 



The main purpose of this Letter is to develop a velocity 
selection theory for LFM, including in the consideration 
anisotropic surface tension effects. We note that this is 
not a just routine extension of the existing theory because 
the diffusion interaction between two interfaces changes 
the problem substantially. 

We discuss the two-dimensional problem of the steady- 
state motion of a thin liquid film during the process of 
isothermal melting of a binary alloy, see Fig.l. We as- 
sume that the diffusion in the solid phases is very slow 
and the concentration c in the liquid film obeys the 
Laplace equation. We introduce the normalized concen- 
tration C = {c — cl) / (cl — cs) with cl and cs being 
the liquidus and solidus concentrations of the equilibrium 
phase diagram at a given temperature. Then the equilib- 
rium concentration and the mass balance conditions at 
the solidification front read 



C = d2K2 



-DdC/dn. 



(1) 



At the melting front the equilibrium concentration is 
changed by the presence of the elastic coherency strain 
energy Q and also the diffusional flux changes because 
in the solid ahead of the melting front the concentration 
is Co which is different from cg: 



C = -6A2 - diKi 



A) = -L>9C/an. (2) 



Here Vn is the normal velocity; D is the diffusion co- 
efficient in the liquid film; K is the curvature assumed 
to be negative for the interfaces in Fig.l; A = (co — 
cs)/{cl — Cs) is the dimensionless driving force; b — 
Y Vl{da / dcY / a? f'l is the dimensionless constant which 
describes the coherency strain energy Q , is the atomic 
volume, Y is the bulk elastic modulus, a is the atomic 
constant, /l(c) is the free energy of the liquid phase per 
atom, f'l is the second derivative of fL{c) at c = cl; di 
are the anisotropic chemical capillary lengths: 

d,{e) = do[l - acos4(0 - e,)], do = ^n/f'l{cL - cs? 

(3) 



with the isotropic part of the surface energy 7 and the 
anisotropy parameter a << 1; is the angle between the 
normal to the interface n and the direction of motion, 9i 
is the direction of the minimum of di{9) for each of the 
interfaces. 

We measure all lengths in the units of the radius of 
curvature of the Ivantsov parabolic solidification front, 
i?2. Introducing a parabolic coordinates system. 



(4) 



we look for a solution of the Laplace equation for the 
concentration field C(?7, [3) in the liquid phase with the 
following boundary conditions: 

c{in = m{P).P)=d2{m2{P), 



C{tj = mW,P) = -bA^-d,{f3)K,{f3). (5) 
The mass balance conditions at the interfaces arc: 

2P(l-A)(7]i(/3)+/3?7'i(/3)) = ~{dC/di^~i^[dC/dl3) (6) 

where the Peclet number is P = VR2/2D and V is the 
steady-state velocity in the y-direction. We present the 
two moving interfaces in the form, 

??2 = l + e2(/3), Vi=VQ + eiW (7) 
and the concentration field in the form, 



(8) 



where Co{ri) = — 6A^(7/ — 1)/ (ryo ^ 1) is the solution with- 
out surface tension (w = ei = £2=0) which satisfies the 
boundary conditions: Co(l) = and 6*0(770) = — 5A^. 
The balance equations at the interfaces, 2P = 2P{1 — 
^)Vo = —dCo/drj lead to the relations for the Peclet 
number and 770: 



P = 



VR2 
2D 



6A2 



1 



2(770-1)' 



Vo^ \ -FT = 



(9) 



These relations can also be obtained in the proper limit 
of Eqs.(27)-(28) of Ref. 01 and give the expressions for 
the radii of curvature for the two interfaces, i?i and i?2, 
for a given velocity V which remains undetermined at 
this stage. 

The surface tension plays a crucial role in the velocity 
selection problem. In order to find the small corrections 
to these solutions due to the small (but singular) surface 
tension effects we should first find the small correction u 
to the diffusional field in the linear approximation with 
respect to e and dK (see, for example, 0|). In this 
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approximation, the field u satisfies the Laplace equation 
and the boundary conditions: 

u{tj = 1, /3) = ^2{I3) = 2Pe2 + d2K2, 



m(?7 = 7?o,/3) = = 2Pei-di/u. (10) 

The balance equations read: 

2P(e2 + pe'^) - ~du/dr]\^^i 



2P{1 - A)(ei + f3e[) = -du/dTj\^^ 



(11) 



The Laplace equation with the boundary conditions, 
Eas. Hl()|l . is easily solved by Fourier transform, with re- 
spect to the variable /3, which we define as 



<?(A) 



1 



/2n 



d/3exp(-iA/?)g(/3) 



dX exp{iX0)g{X). 



Then, 



— — _ ^1 sinh[A(7? - 1)] ^2 sinh[A(7? - tjq)] 

> sinh[A(7yo - 1)] sinh[A(77o - 1)] ' ^ ^ 

Finally, Eas. (fTT|l read 

r,pde^ ^ IT WTcosh[A(??o - 1)] 

dX sinh[A(77o — 1)] sinh[A(77o — 1)] 



2P(^ - A)^ ^ ^icosh[A(?7o~l)] _ ^2 

^ ' dX sinh[A(r/o - 1)] sinh[A(7?o - 1)] ' 



(13) 



Eliminating ei from the first of these two equations, 
2Per = 2P(ei)' sinh[A(?7o - 1)] + 2FeJcosh[A(77o - 1)]- 



diKi + d2K2 cosh[A(77o - 1)], 
we find from the second one. 



(14) 



2P(e^)" - 2Pei ^ ^2/^2 



{diKi)' 



sinh[A(77o - 1)] 



{d2K2) 



,cosh[A(7;o - 1)] 



sinh[A(77o - 1)] 



(15) 



Here ' denotes derivatives with respect to A. 

These equations are convenient for the further anal- 
ysis. We are interested in the motion of the thin film, 
which means that (770 — 1) ~ A << 1. In this case we 
can expand the hyperbolic functions for small values of 



their argument almost everywhere. This leads to ei w £2, 
Ki w K2 and then, from Ea. (|15() we find in direct space: 



2PA 



■-PK2 



(16) 



A few remarks are in order. First of all, this is somewhat 
the expected result because for the thin films the result- 
ing equation should be local and, in principal, can be 
derived in the framework of boundary layer techniques. 
However, the used approximation for small arguments of 
hyperbolic functions breaks down in the small vicinity (of 
the order of A) near the singular points in the complex 
plane, (3 = ±i. We will discuss this point later. Finally, 
we can return to Cartesian coordinates because the ex- 
pression for the curvature is much simpler in this repre- 
sentation. In the linear approximation (l + /3^)e2 ~ C,2(x) 
where C2 {x) is the correction to the parabolic solution in 
the Cartesian coordinates. Then Ea. H16|) reads: 



d<2 
dx 



di + d2 
2PA 



(17) 



One finds the first regular corrections to the parabolic 
shape by replacing di — d2 ^ dp and K2 = — l/[i?2(l + 
2;2-j3/2j -j-Qj. ^j^p parabola: 



C2(a;)-a/[A(l + x2)i/2], A « 1, 



(18) 



where a = do/{PR2) is the usual stability parameter 
which appears in such problems. 

In the vicinity of the point x = —i where the curva- 
ture becomes singular and anisotropy, Eq.Q, becomes 
important we rescale the variables 



0, C2 



aF, 



z + dF/dz, (19) 



with a small anisotropy parameter, a << 1. Then 
Eq.((T7I) reads 



dF 

dz 



Aa5/4 



1 + d^F/dz^ 
(2t)3/2 



(20) 



where the anisotropy factor is A{t) = 1 — [exp(— 4i6'i) + 
exp(— 4i6'2)]/T2. This equation is similar to the limit of 
the kinetic dendrite in Eq.(9.1) of Ref . [l(l| . This simi- 
larity is purely formal because the real interface kinetic 
effects are not included in our description. Instead, the 
strong diffusion interaction leads to such a selection prob- 
lem. The same arguments as in [lo| give the selection 
condition 



A << 



(21) 



Moreover, the selected growth direction lies in between 
the directions of the minima of c?i {9) and ^2 {0) for the two 
interfaces, 9i = —O-i, because the coefficient appearing in 
A{t) must be real jl(l |. 

While the regular correction given by Ea. H18|l is valid 
for A << 1, the restriction for the validity of inner 
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Ea. (|20|l and for the scaling relation, Ea. (|21() . is much 
stronger. A << ^Ja « 1. In the opposite limit, 
A >> -y/a, we can derive another local inner equation 
where the important scale |x + i| ~ \fot is much smaller 
then A. The argument of the hyperbolic functions be- 
comes large and we can neglect the term K\ in Ea. (|15|l . 
Then, in the vicinity of the singular point x = —i in the 
complex plane, this equation reads: 

C2ix)^-d2K2/P, (22) 

which is the usual inner equation for a single dendrite in 
the one-sided model With the same rescaling as in 
Eq.ldni), we find 

~ aV^ ^ ' (2t)3/2 ' ^ ' 

where the anisotropy factor is A{t) = 1 — 
2exp(— 4i6'2)/T2. This corresponds to the selection con- 
dition 

a = ^7*-Q^^/^ A>>yS (24) 

and the selected growth direction corresponds to the min- 
imum of d2{6), 02 = 0. Whereas Ea. (|23|) and the selec- 
tion condition Eq. (|24|) are valid not only for small A 
[lo| , the regular correction in the form of Eq. (|18|l requires 

A«l. 

The selection conditions, Eas. (|21l24|l . together with 
the relations for the Peclet numbers, Eq.®, solve the 
posed problem of pattern formation for two combined 
nearly parabolic fronts: 

9 D 

V=—P^{A)a* (25) 

"0 

and 

R2^do/[P{A)a*], Ri=R2/{l-Af. (26) 

While these results are formally similar to the free den- 
dritic case, the selected stability parameter a* scales as 
a^/^ only for A >> -^a. Otherwise, the other scaling 
relation, Ea. pH) . holds. In principal, it could be con- 
ceivable that the strong diffusional interaction between 
the front leads to the selection even without anisotropy. 
However, our analysis of the inner equation does not sup- 
port this hypothesis. The Peclet number is P{A) ~ bA 
which reflects the fact that the coherency strain energy 
plays a crucial role in the LFM mechanism. This pa- 
rameter b is usually small but the melting process by the 
LFM mechanism is controlled by the fast diffusion in the 
liquid, whereas the the process with only one melting 



front would be controlled by the very slow diffusion in 
the mother solid phase. As we already noted, with the 
help of the LFM mechanism the system chooses a more 
efficient kinetic path to relax to the equilibrium state. 

Finally, we estimate the velocity V from Ea. (|25|) and 
the thickness of the film R2A from Ea. H2(j|) using char- 
acteristic values of the parameters: D ^ 10~^cm2/s, 
do ~ lO-'^cm, b - 0.05, A ~ 0.05 and a* - lO'^. This 
leads to ^ lO^^cm/s and R2A ^ lO^^'cm, which qual- 
itatively agree with typical values in LFM experiments. 

In conclusion, we developed a selection theory for the 
process of liquid film migration where the strong diffu- 
sion interaction between melting and solidification fronts 
plays a crucial role. This process is very important in 
practical applications, in particular during sintering in 
the presence of the liquid phase Q,!^. However, despite 
its practical importance, experimental investigations of 
this process are far from the level of accuracy of the model 
experiments in classical dendritic growth. Our approach 
extends selection theory developed for dendritic growth 
to LFM and we hope that our results will stimulate fur- 
ther theoretical and experimental investigations in this 
very interesting field. From the theoretical side it would 
be a challenge to attack this problem by a direct numer- 
ical approach, for example, by means of the phase-field 
model. 
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